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A measurement of the branching ratios of D + and hadronic decays to four-body 

final states containing a Ks 
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We have studied hadronic four-body decays of D + and Df mesons with a Ks in the final 
state using data recorded during the 1996-1997 fixed-target run at Fermilab high energy pho- 
toproduction experiment FOCUS. We report a new branching ratio measurement of F(D + — > 
K s K--k + -k + )/T{D + -» K s n + ir + TT~) = 0.0768 ± 0.0041 ± 0.0032. We make the first observation 
of three new decay modes with branching ratios F(D + — ► KsK + -k + -k~)/F(D + — * Kstv + -k + tt~) — 
0.0562±0.0039±0.0040, T{D + -» K S K + K--k + )/Y{D + -> Ksn+TY+n') = 0.0077±0.0015±0.0009, 
and T(Df -» K s K + ir + iT-)/r(Df -» K S K~ -k + -k+) = 0.586 ± 0.052 ± 0.043, where in each case the 
first error is statistical and the second error is systematic. 

PACS numbers:13.25.Ft, 14.40. Lb 



To understand hadronic decays of a heavy quark sys- 
tem one needs to address final state interactions. These 
become more complicated as the number of final state 
hadrons increase, since the signatures of the weak de- 
cay of the heavy quark are masked by the hadronic de- 
grees of freedom. Theoretical predictions are still limited 
to two-body decays, which have been analyzed exten- 
sively in the theoretical literature For example, 
Bauer, Stech and Wirbel fl have used a factorization 



approach. Bedaque, Das, and Mathur [El and Kamal, 
Verma and Sinha have used heavy quark effective 
theory. Much less is known about four-body hadronic 
decays than about two or three-body decays. More ex- 
perimental data on higher multiplicity modes is essential 
to improving our understanding of the decay process in 
heavy quark systems. In this letter we report on D + and 
branching ratios into four-body final states involving 
a Kg- We measure the D + decay rates into K$K ~n + ir + , 
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KsK + it + it~ , and K$K + K~tt + relative to Kstt + tt + tt" 
and the decay rate of Df — > KsK + ir + ir~ relative to 
Df — > i^5i ; L _ 7r + 7r + (throughout this letter the charge 
conjugate state is implied). Among these final states 
only the KsK~tt + it + final state has been observed pre- 
viously Q. These final states have several interesting 
resonant contributions which will be the subject of a fu- 
ture report. 

We collected the data for this study during the 1996- 
1997 fixed-target run of the photoproduction experiment, 
FOCUS at Fermilab. This experiment utilized a for- 
ward multi-particle spectrometer to study charmed par- 
ticles produced by interactions of high energy photons, 
(-Ely) w 180 GeV, with a segmented BeO target. Charged 
particles were traced by silicon microstrip vertex detec- 
tors. These detectors provide excellent separation be- 
tween the reconstructed production and decay vertices. 
The vertex resolution is approximately 6 /jm in the trans- 
verse direction and 300 /im in the longitudinal direction. 
Three multi-cell threshold Cerenkov detectors were used 
to identify charged hadrons. 

We reconstructed the K$ candidates using the decay 
Kg — > 7r + 7r~ and calculated the error on the K$ mass 
for each candidate. We required the Ks mass be within 
3cr of the nominal Ks mass which rejects nearly all of the 
e + e~ pair background. To reduce the A — > pir~ back- 
ground we applied a Cerenkov particle identification cut 
of 5 units on the difference in the log likelihoods between 
the pion and proton hypothesis || on the Ks daughter 
track with higher momentum. 

We selected the final states using a candidate driven 
vertexing algorithm || . A secondary vertex was formed 
from the reconstructed tracks and the momentum vec- 
tor of the charm candidate was used as a seed track in 
finding the primary vertex in the event. We required 
that the primary and secondary vertices be formed with 
a confidence level greater than 1%. The significance of 
separation between the primary and secondary vertex is 
called £/&£. We required l/ai > 9 for D + candidates and 
I jut > 7 for £>+ candidates. We also required the proper 
decay time be less than 5 times the candidate particle's 
lifetime [ |l0[ . The vertexing algorithm provides two esti- 
mators of the relative isolation of the vertices. The pri- 
mary vertex isolation (ISOl) estimator is the confidence 
level of the hypothesis that a track in the secondary ver- 
tex was also in the primary vertex; the secondary ver- 
tex isolation (IS02) is the confidence level that another 
track originates from the secondary vertex. We required 
ISOl be less than 1% and IS02 less than 0.1%. For the 
D + — > KsK~tt + tt + decay mode, this removes 80% of 
the background and retains 70% of the signal. Further, 
we required that the decay vertices be out of target ma- 
terial, eliminating backgrounds from interactions which 
are induced by particles from the primary interaction or 
from conversions of spurious photons. 

In addition to the vertexing requirements, we used 



Cerenkov particle identification cuts. For the charged 
kaon candidates we required the kaon hypothesis be fa- 
vored over the pion hypothesis by 2 units of likelihood 
(for the D + — > KsK + K~ir + , we required 2 units for the 
faster kaon but only 1 unit for the slower kaon) . 

Finally, the charm candidates were required to have a 
momentum greater than 30 GeV/c, which removed com- 
binatoric non-charm backgrounds. For all decay modes 
we selected cuts by maximizing the figure of merit defined 
as A/' s 2 ignal /(A/' s ignai + A"background), as well as minimizing 
reflection background. 

We investigated several possibilities for contamina- 
tion of either signal region due to reflections from 
partly reconstructed charmed particles. We find no ev- 
idence in the data for an enhancement due to A+ — > 
Kspir + ir~ where the p is misidentified as a if or due 
to A+ — » A°7r + 7T + 7r _ where the A is misidentified 
Ks- We also investigated contamination of the 
Df — > KsK~ir + ir + and — > KsK + tt + it~ signal re- 
gion due to Tt/K misidentification of the decay D + — > 
AT,j7r + 7r + 7r~. This reflection is incorporated in the fit 
(see Figs, [l] and ||). Double misidentification between the 
KsK~ir + Tr + and KsK + t: + it~ decay modes was found to 
be negligible. 

In order to minimize systematic effects we chose nor- 
malization channels which have similar topologies to 
the signal modes, and we applied the same vertex re- 
quirements. Fig. ^ shows the invariant mass distri- 
butions using the D + selection cuts and Fig. shows 
the invariant mass distributions using the Dj selec- 
tion cuts. We parameterized the signal with a Gaus- 
sian and the background with a quadratic polynomial 
plus a reflection shape. For the D + — > Kstt + tt + 7t^ 
normalization mode we used a linear background. The 
reflections in Figs. |(a), 0(b) and Figs. |(a), |(b) 
from D + — > Ks7r + ir + TT~ and Fig. 0(c) from D + — > 
KsK~Tr + Tr + (KsK + ir + ir~) occur due to the misidenti- 
fication of a 7r as a K. The shapes were determined by 
Monte Carlo simulations and the reflection amplitudes 
are free parameters in the fit. 

The final states in this study have resonant substruc- 
ture in their decay modes which may affect the recon- 
struction efficiency. We generated a Monte Carlo for 
each final state using an incoherent mix of various sub- 
resonant decay modes and a non-resonant contribution. 
For the decay D + — > Ks^ + tt + tt^ we used its known sub- 
decay modes to determine the efficiency. Mini-Monte 
Carlo studies showed this efficiency is consistent with 
that of allowing interferences between the states. For the 
other decay modes, where little is known about the res- 
onant substructure, we determined the efficiencies based 
on which sub-decay we found to be dominant. The de- 
cay mode used for the KsK~tt + tt + (KsK + K~ir + ) fi- 

nal state was K K* + {Ks<j)-K + ). The K s K + tt + it~ fi- 
nal state has contributions from many sub-decay modes; 
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FIG. 1. The invariant mass distributions with D + selection 
cuts, (a) shows the final state KsK~tt + tv + , (b) KsK + 7r + 7r~ , 
(c) KsK + K~n + and (d) Kstv + tv + tt^ . The data are indicated 
by points with error bars and the solid line is the fit. The 
shaded regions are charm reflection backgrounds. 



therefore we assumed a non-resonant substructure to ob- 
tain the efficiency for this final state. With these Monte 
Carlo efficiencies, we measured the relative branching ra- 
tios which are summarized in Table |. For each of the 
signal modes we also obtained efficiencies with different 
resonant and non-resonant substructure combinations. 
We used the spread of the various resonant and non- 
resonant efficiencies to determine the systematic error on 
the efficiency calculation for each state. We estimated 
the relative uncertainties in efficiencies due to uncertain- 
ties in the branching fractions into resonant sub-states 
to be 3.9% for D+ -> K S K -7T+7T+, 6.7% for D+ -> 
KsK+ir+ir-, 9.8% for D+ -> K s K + K-<k + , 4.6% for 
Df -> K s K-it + it + , and 4.6% for Df -> K s K + tt+^- . 

A study of the stability and behavior for each branch- 
ing ratio was performed using variations of our analysis 
cuts; we found no bias from the choice of analysis cuts. 
Further, we split our data into independent subsamples 
based on D momentum and the different run periods in 
which the data were accumulated. This technique is de- 
scribed in detail in reference jll]] . We found no systematic 
uncertainties from splitting our data. 

The systematic uncertainty for each branching ratio 
includes efficiency dependencies from sub-resonant states 
and from variation of the fit parameters (mostly due to 
background parameterization variation). 

To conclude, we have improved the previous mea- 
surements of the D + — > KsK~tt + it + branching ratio 
and have made the first observation of the D + (Df) —> 
KsK + tt + tt~ and D + — > KsK + K~n + decay processes 
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FIG. 2. The invariant mass distribution with Df selec- 
tion cuts, (a) shows the final state KsK~-k + -k + and (b) 
KsK + tt + -k~ . The data are indicated by points with error 
bars and the solid line is the fit. The shaded regions are 
charm reflection backgrounds. 





Decay Mode 


A/signal 




D+ 


— > KsK~Tr + Tr + 


670 ±35 


0.0768±0.0041±0.0032 


D+ 




469 ±32 


0.0562±0.0039±0.0040 


D+ 


-> K s K+K-ir+ 


35 ±7 


0.0077±0.0015±0.0009 


D+ 


— » Ag7T + 7T + 7r~ 


11590±121 1 


Df 


— » KsK + n + Tr~ 


476 ±36 


0.586±0.052±0.043 


Df 


— > KsK~TT + Tr + 


837 ±38 


1 



TABLE I. Measured branching ratios. r ro i is the branching 
ratio relative to D + — > Ksir + n + Tv~ for the D + modes and 
Df —> KsK~ 7r + 7r + for the Df modes. The errors on the 
branching ratios are statistical and systematic, respectively. 



using a high statistics sample of photoproduced charmed 
particles from the FOCUS experiment at Fermilab. 
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